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The spectral function of one hole in different magnetic states of the one-dimensional t-J model 
including three-site term and frustration J' is studied. In the strong coupling limit J — > (corre- 
sponding to U —* oo of the Hubbard-model) a set of eigenoperators of the Liouvillian is found which 
allows to derive an exact expression for the one-particle Green's function that is also applicable at 
finite temperature and in an arbitrary magnetic state. The spinon dispersion of the pure t-J model 
with the ground-state of the Heisenberg model can be obtained by treating the corrections due to 
a small exchange term by means of the projection method. The spectral function for the special 
frustration J' — J/2 with the Majumdar-Ghosh wave function is discussed in detail. Besides the 
projection method, a variational ansatz with the set of eigenoperators of the t-teim is used. We 
find a symmetric spinon dispersion around the momentum k — 7r/(2a) and a strong damping of the 
holon branch. Below the continuum a bound state is obtained with finite spectral weight and a very 
small separation from the continuum. Furthermore, the spectral function of the ideal paramagnetic 
case at a temperature ksT 3> J is discussed. 
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I. INTRODUCTION 



> 

■ The understanding of spin-charge separation is a central point in the physics of low-dimensional electronic systems. 
That issue is most clearly seen in one dimension where the exact solution of the ID Hubbard modelu reveals that the 
low-energy physics is dominated by decoupled, collective charge and spin excitations (also called holon and spinon, 
respectively). The idea that the spin and charge degrees of freedom separate has also been proposed to explain the 
properties of 2D cuprate superconductors .□ However, even in one dimension, the Bethe ansatz solution did not vest 
give a complete answer for the spectral function A(k, u>). Only for U — ► oo an exact expression for A(k, lu) was foundEl 
due to the factorization property of the wave function accordinahto Ogata and Shiba.cl The results for the insulating 
half-filled caseQ could also be generalized to other filling factors.B 

The spin-charge separation was observed in ARPES measurements of the one dimensional, dielectric cuprate 
, SrCu02.Q The spinon and holon branch of the spectral function were seen, in contrast to the analogous experi- 
ment of one hole in the Cu02 planeQ where spin and charge are coupled and the spin polaron quasiparticle has a 
dispersion proportional to J as proposed theoretically in Refs. [8|-|l0| The ARPES spectra in SrGu02 were analyzed 
using the pure t-J model. On the other hand, many ID compounds, like for instance CuGe03,U are characterized 
by frustration in the magnetic subsystem which may lead to a gap in the spin excitation spectrum. Fox,the special 
frustration J'/ J = 0.5 and in the limit J — > an analytic expression for A(k,u>) was derived recentlyji3 under the 
assumption that the wave function factorizes. Besides the frustration, also temperature effects are important as it 
was observed in ARPES measurements on Nao.96V2O5.t2l So, there is a clear need to study the spectral function 
systematically under the influence of frustration and temperature and to derive analytic expressions. 

The present work focuses mainly on the effect of frustration and temperature _pn the spectral function in the 
insulating case. For that, we rederive first the exact solution of Sorella and Parolatj in a straightforward way using 
Green's function technique. That is possible due to our finding of a set of eigenoperators of the Liouvillian (Sec. Ill) 
in the strong coupling limit J — ► 0. As a consequence, our derivation is applicable for any mag netic state and any 
temperature in that limit. Especially, one can show that the result which was derived in Ref. [12|does not depend on 
the assumption that the wave function factorizes. We present analytic expressions for the spectral function of one hole 
in several magnetic states: (i) the ground-state of the antiferromagnetic Heisenberg model, (ii) the Majumdar-Ghosh 
wave functionllil at the special frustration J'/ J = 0.5, and (iii) the ideal paramagnetic state at temperatures much 
larger than the exchange energy fc^T ^> J. 



For large finite coupling we compare two methods to account for corrections J cx (t 2 /U), namely the projection 
method and a variational ansatz using the set of eigenoperators of the i-term. We show that the former method yields 
a reasonable description of the spinon dispersion in the pure t-J model (Sec. IV) and an approximate result for the 
spectral function of the Majumdar- Ghosh model. For J' = J/2 it misses the bound state below the continuum which 
is obtained by the more accurate variational method (Sec. VB). The bound state has a finite spectral weight but a 
very small separation from the continuum. Both methods show in the Majumdar- Ghosh case that the low energy 
region for momenta k between tt/2 and it (lattice constant a — 1) will be filled with states, that the spinon dispersion 
(i.e. that collective excitation corresponding to the lower edge of the continuum) becomes symmetric around tt/2, 
and they indicate an overdamped holon branch. The damping of the holon branch is extremely large for very high 
temperatures (Sec. VI). 

Before presenting our results let us shortly discuss the different understandings of the term "spin-charge separation" 
as it can be met in the literature. The naive picture means that the low energy effective Hamiltonian may be written 
as 

H = H h + H s , [H h ,H s }=0, (1) 

and the electron operator is the product 

= s i(T hl , (2) 



c 



where spinon s and holon h can be basically regarded as free particles. Then the normalized (i.e. J A(k,uj)duj = 1) 
spectral function is 

A(k, ^) = \Y. 2 f^ 6 l u - e "( fc " 0) " e ^ ' ( 3 ) 
Q 

where f(Q) = #(§ — \Q\) is the Fermi distribution function of spinons, 9{x) is the Heaviside step function, L is the 
number of sites, e^, e s being holon and spinon energies, respectively. _ . 

However, the naive understanding is not that one which is realized in ID electron systems.E3 There, it was found 
that the eigenstates factorize in the limit U — > oo in the form 

tp(xx, ...,x N ,yx. ..y M ) = *Psf(xi, ■ ■ -,x N )(j) H {yi . ..y M ), (4) 

where xi, . . . ,xn are the spatial coordinates of the N electrons on a L-site ring, and the y\ . . .yu 'coordinates' 
label the position of the spin-up electrons on the squeezed Heisenberg ring, i.e. on the N occupied sites. The ipsF 
is a spinless fermion state, and 4>h is an eigenstate of an iV-site Heisenberg Hamiltonian with periodic boundary 
conditions. The product form of equation (0) should not be interpreted as a trivial decoupling between charge and 
spin. In fact, the momentum of the spin wave function imposes a twisted boundary condition on the spinless fermion 
wave function. As a result, the Fermi distribution function f(Q) in (||) will be replaced by a function Z(Q) that is the 
expectation value of a chain of spin operators that has to be determined from the pure spin systemu (for the details see 
Sec. III). The singularity of Z(Q) produces additional peaks in Alk,uj). This correct answer for the spectral function 
may be understood as a manifestation of the phase string effectJij It means that spinon and holon interact with each 
other via a nonlocal phase-string. Instead of (0) we should write 



= s ia h\ exp 



l>i l>i 



It should be noted that the phenomenon of spin-charge separation is not restricted to the limit U —*rPO m the ID 
Hubbard model. At any finite U the spin and charge fluctuations propagate with different velocities!] That means 
that after some time the spin and charge degrees of freedom will be separated in space. But there is no analytic 
solution for the spectral function of the Hubbard model at arbitrary values of U and also the present calculation 
treats terms of order (t 2 /U) as a perturbation. In that sense we will understand here spin-charge separation as a 
manifestation of the factorization property (^) in the spectral density. Sharp maxima in the continuum correspond 
to collective excitations whereas possible bound states indicate special eigenfunctions with a strong coupling between 
spin and charge. Another possible effect of additional terms in the Hamiltonian is the broadening (i.e. the damping) 
of the collective excitations. 
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II. MODEL AND SPECTRAL DENSITY 



To describe the low energy physics of compounds with a ID electronic structure it is sufficient in most cases to 
take into account only that band which is closest to the Fermi energy (see for instance Ref. [n]). Treating the on-site 
Coulomb interaction explicitly, one obtains the well known ID Hubbard model. In the present calculation we restrict 
ourselves to the strong coupling limit U ^> t where we may project out the subspace of doubly occupied sites, and for 
the lower Hubbard band we obtain the effective Hamiltonian 

H = i+j + t 3l (5) 

where 

t = -tY J XfX^ g , (6) 



J \ - Y af3 Y l3a 

2 2^ A * A *+i ' 



Xr P X^ , (7) 



I _. \ " yaO y^a Y 0/3 r Q \ 
*3 - *3 2^ A i X i+g X t+2g i (°) 

and a, (3 =TI; 3 are the nearest neighbors g = ±1. The Hamiltonian is valid near half filling {X^ + + X^~ = 1). The 
parameters J and £3 are connected with the original values of the Hubbard model by 

J = 4t 3 = 4t 2 /U, (9) 

but the Hamiltonian (^) is more general, if we relax the condition (^). It may be derived directly from the more 
realistic three-band Hubbard model. Then, the i3-term often becomes negligible and one obtains the t- J model. The 
Hamiltonian (|5|) is written in terms of Hubbard projection operators that act in the subspace of on-site states 

Xf= \a,i)(0,i\ , a,/? = 0,U,2. (10) 

They are related with bare fermionic and spin operators through 

X? = cf (1 - rH,-a) , Xf = -aa^m^ , (11) 



x+- = s+ = c J tCi , , xr = 



c*, T c a , X- = \ + \ (ct T c i>T - c\ d c iA ) = i + aSt , (12) 
with a = ±1. Other relations are easy to obtain with the use of the main property of Hubbard operator algebra 

Xfxf = 5 [3l X^ , (13) 
that follows immediately from the definition (|l0|). The commutation relations for operators on different sites are 



fermionic for operators that change the number of particles by odd integers, like (11), and bosonic for others. In the 
presence of frustration in the magnetic system, which is discussed for instance for CuGe03, the t-J Hamiltonian may 
be generalized by inclusion of the J'-term 

i' = yE x ? x & ■ ( 14 ) 

Our aim is to calculate the one-particle two-time retarded Green's function G(k, to) and the spectral density of one 
hole in the magnetic state 

A(k,u>) = ImG(k,uj + iO + ), (15) 

that is roughly proportional to the ARPES signal intensity. We define 
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2nS(k- k')G{k,u) = {{X?\Xl?)) ee -i J dte^-VaxgOfoXptf)}) , (16) 

where 

+00 



({X?,Xp}) = 2*8(k-k'), 

and where {...,...} means the anticommutator. The expectation value denotes the thermal average over a grand 
canonical ensemble: 

(...) = Q- x Sp [e-^~^)...], Q = Sp e -MB-i>N). (17) 

Here Sp implies taking the trace of an operator, N is the particle number operator, j3 — (fcT)" 1 is an inverse 
temperature, and \i represents the chemical potential. The time dependence of the operator B(t) is given by B(t) = 

e it{H-iiN) Be -it(H-/j.N) _ 

III. EIGENOPERATOR AND HOLON DISPERSION 

Let us consider first the limit U — > 00 in the Hubbard model or J, £3 — > in (|5|). Then only the i-term 

t = -tY,xfx^ g , 

is nonzero. Note that it is a true many-body Hamiltonian due to the constraint of no double occupancy, as we see 



from Eq. (11). We introduce the set of operators 

Y X~iX%™...X%£°;X%° r , ff = sign(r), (18) 



Qi , . ■■■a r 



for which 

[Vm.r, t\=t (U m ,r— 1 + V m ,r+l) (19) 



holds at half filling. Any operator (|18|) can be considered as a string operator of a certain length consisting of a hole 
and an attached string of spin flips. Such a string can be produced by creating a hole in the Neel state and applying 
several times the kinetic energy (^|) which jcijeates misaligned spins. Similar string operators were used to describe the 
spin polaron quasiparticle in the 2D case.E3o We make double Fourier transform 

v k , q = V2 J2 e~ ikm ~ iqr v m ,r , (20) 

m,r— — oo 

and see that 

[wfc, g ,t] = 2tv ktq cosq . (21) 

The interpretation in terms of the string operator ([^J) is quite easy. We see that only the right end of the operator 
( |l8| ) was influenced by the t-term. Therefore, we may identify the right end of v m ^ r with the holon excitation. In the 
next Section, it will become clear that the left end of v m . r may be connected with the spinon. 

The operators Vk, q ( p0| ) are eigenoperators of the Liouvillian L of the problem, where CA = [H, A]. Note that it is 
one of the rarest, if not the unique case in many-body physics that the explicit form for a set of eigenoperators can 
be given. From ( |2l| ) we see that the equation of motion for the corresponding string operator Green's function closes 
and it has a simple pole form 

« v *«\ v «.*)) = ^t-a^s'P ' = 87r ^ k k 'M q q '^ k + ( 22 ) 
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where the spectral weight is the expectation value 



Z{q + *) = \ £ e-^(O r ) , (23) 



of a chain of AT-operators 



Q \ * y(Tai yaift2 yOir—lOir ya r (7 

r ~ / j ^m+g ■ ■ ■ ^m+r-g^m+r 



ai,...,a r i<7 



— (2S m S m +g + — )(2S)T[+gS m +2g + — ) • • ■ (2S m +) — g S m + r + — ) , (24) 

to be calculated for the pure spin-system without any hole. The expectation value (fir) of (|2i|) cannot depend on 
the starting point m due to the translational symmetry of the problem. The operator ( p4[ ) was introduced in Ref. ^ 
and explicit values on a 26-site Heisenberg ring were given for T = 0, when (. . .) becomes the average over the 
ground-state. Asymptotically, the following behavior was found 

(flj) -> -^Re \Ae tnl/2 } (25) 

which leads to a square root singularity of Z{Q). Using additionally the exact values (Oo) = 1 and (Ox) = 1 — 21n2 
the following formula may be derivedH 

Z{Q) = (-0.393 + 0.835/ VcoIq) 0(| - \Q\) ■ (26) 

For the hole Green's function ( |l6| ) we have 

J_ v 2ttJ_ 7[ 2tt k > q 7T uj-2tcosq 

and the spectral density is obtained in the way 

r+ n dO 

A(k,w)= —Z(Q)5[u) + 2tcos(k-Q)]. (28) 

J-7T T 

That gives the exact answer in the strong coupling limit (J — > 0) where only the holon dispersion e^(q) = — 2tcosq 
is important. But the ratio J' / J may be arbitrary and ( |28|) is not only exact in the Heisenberg case with Z(Q) from 
(p6|). Instead, from our derivation follows its validity for arbitrary magnetic states and it is not restricted to zero 
temperature. Then, however, Z(Q) is different. In the following we will give exact results for i) the Majumdar- Ghosh 
wave function at the special frustration J'/ J = 0.5, and ii) the ideal paramagnetic case at fc_eT 3> J. Two cases 
are quite trivial, namely the saturated ferromagnetic case and the classical Neel state. The former one leads to 
Z(Q + 7r) ql.5(Q) and a spectral function like for free fermions, whereas the latter case leads to the Brinkmann-Rice 
continuumEl (Z(Q) = 1/2). A magnetic state inbetween the classical Neel and the Heisenberg case could, in principle, 
also be considered, for which the one-hole spectral function was derived in Ref. ^ treating the spin-fluctuations as 
perturbation. 



IV. SPINON DISPERSION 



In real systems the ratio J ft is roughly 0.3. It means that they are in the regime of strong coupling and the above 
consideration correctly describe the largest energy scale oc t. Now, we want to estimate the corrections that arise 
from other terms of the Hamiltonian (Eh. First we note that Vk, q are eigenoperators for the t3-term 

[vk,q, h] = -2i 3 Wfc, ? cos 2q , (29) 

which leads to the replacement thi,q) — * £h(q) + 2t% cos2<7 in the denominator of (p2|). The commutation with J gives 
(see Appendix A for the details) 
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J. 



V k,q, J = - JcOs(fc - q)Vk,q + ^K,g + Wfc'g) , 



(30) 



where u^. and i^' are Fourier transforms of the operators 



^m,r / j m rn — g m+g 



' -A-m+r-l-A- 



m+r i 



(31) 



E 

7,Ql ,. . . 



• m-\-r — g m-\-r-\-g rn-^-r * 



(32) 



It is impressive that terms, which come from the commutation of "inner" X n operators in v m ,r with n between the 
points to and m + r cancel each other and only the terms coming from the ends remain. The term v'^ r presents 
a distortion of the right end of v m , r by means of the exchange part and may be interpreted as the loss of magnetic 
energy due to the presence of a holon. On the other hand, the term v' m r , with a distorted left end will be shown to 
give rise to the spinon dispersion. We really observe the "separate" motion of the holon that is represented by the 
right end of u m ,r and of the spinon that is the left end of u m>r . The holon motion is governed by the i-term and the 
spinon motion by the J-term. We put the word "separate" in quotes because the motion remains correlated due to 
the set of "inner" X n operators, connecting the ends of v m _ r . 

We need an approximate approach to account for v' k + v'l q . For this purpose we use the projection technique 



J k,q 



J k,q 



{v' k , q +v'l,qA,q}) 



Vk,q ■ 



({«*,«> 4,,}) 

Now, the Green's function for the string operator has the form 

t _ 87T 2 <5(fc - k')5{q - q')Z{k -q + Tl) 

\\ v k,q\ v k>,q'// - UJ -2tcosq + 2t 3 cos2q-e s (k-q)' 

where e s (k — q) is defined by the equation 

({ [«fc, 9 . J] . 4',q'}) = - k')S(q - q')e s {k ~ q)Z{k - q + tt) . 

The contribution of the v'^ term to the spinon dispersion is determined by an expression of the form 



(33) 



2 ^m-\-r,m , -\-r f (^r.r'} ' 



(34) 



(35) 



(36) 



where the precise order of the X-operators in r2" r / can be easily inferred from (|32|) and is given in the Appendix 
A. There, it is also shown that for slowly decaying spin correlation functions (as in the present case, see (p5h) the 
correlation functions (f2" r ,) can be approximated to be a function of r — r' only, in the way: 



(37) 



That leads to a constant shift of the energy e s as the only effect of v'^ r which will be neglected further on. The contri- 
bution of the v' m r term can be written analogously to (pq), defining the spin correlation functions (ft' r r ,) — (£l' r _ r i ). 
The correlation functions (f2J ) differ from (Sli+i) only by the exchange of two A-operators. Therefore, for large I, 
we may expect that 



1+1/ 



(38) 



That leads after Fourier transformation to the contribution of the v' m r term to the spinon dispersion (Appendix A). 
Together with the contribution of Vk, q we obtain for the spinon dispersion 



e s (Q-Tc)Z(Q) 



cos Q 



Z(Q) 



:Ko>-<ni> 



sin Q 



(Ik 
2^ 



Z(k) ( cot 



Q 



cot 



Q 



, (39) 



and the hole spectral function becomes 



G 



A(k,w)= [ — Z(Q)8 [u + 2tcos{k- Q)-2t 3 cos 2(k-Q)-e s (Q ->*)]. (40) 

J-k 7T 



The curve that we obtained for e s with the formula (|39|) is close to 

e s (Q - Ti") « uJcosQ , aw 2. (41) 

as shown in Fig. l._The functional form ( pd| ) is consistent with Bethe-ansatzi and field-theoretical considerations.lil 
(Sorella and ParolaEI derived a contribution Jn/2cosQ w 1.6JcosQ.) We tested it also by comparing the first two 
terms of Fourier expansion of the product cosQZ(Q) with (Sl' ) and (Sl[ n) that give values for a in (^Tj) of 2.1 or 
1.8, respectively (for the pair correlation functions we took the data of Ref. [|). Therefore, we are using the simplified 
formula ( f4l| ) instead of (^) in the following analysis of the spectral density. We have checked that the differences are 
negligible. 

The spectral density (for £3 = 0) is shown in Fig. 2. One can clearly distinguish between the spinon and holon 
features at the lower edge of the spectral density dispersing at an energy scale oc J (from fc = 0tofc = 7r/2) or ex t 
(from k = tt/2 to k =r-f)- At k = k*, which is determined by tcosk* = J, another holon branch splits off the lower 
edge of the spectrumEj and disperses towards k — at an energy scale of t (and a corresponding holon branch splits 
off the upper edge of the spectrum). For k values inbetween and k* one has three peaks in the spectral function 
(one spinon and lower and upper holon branch). One can easily imagine the situation in the doped case. Then the 
spinon and holon branches start at the Fermi energy with two different velocities.Eil 

In contrast to the 2D caseja there is no separate bound state at the lower edge of the spectrum indicating that there, 
are only collective spin and charge excitations. Most of those features were also observed in the ARPES experiment .til 
In the naive picture of spin-charge separation (^) the spectral density would have square root singularities only either 
at the lower or at the upper edge of the spectrum. In Fig. 2, however, there are additional holon branches due to the 
square root singularity in Z{Q) (see also Ref. |l5|). For J — * the spinon feature in the spectral density, i.e. the lower 
edge of the spectrum, becomes completely flat between and 7r/2. The corresponding pictures were already given in 
Ref. |. Fig. 2 agrees also qualitatively with the finite cluster results.B 



V. MAJUMDAR-GHOSH MODEL 



We have shown that our approach is applicable for any magnetic state for J — > 0. Now, we are going to present the 
spectral function of one hole in the t-J-J' model with the special frustration J' = J/2 (called here Majumdar-Ghosh 
(MG) model for simplicity). In that case rigorous analytic results may be obtained since the ground-state wave 
function of the MG spin Hamiltonian is exactly known.Ej It is the combination of two simple dimer states 



*MG = (*i + $2)/ V2 , (42) 

where 

+00 +00 
*i= JJ [2n,2n+l] , $2 = JJ i^-l,2n], 

n — — 00 n=— 00 

and the singlet bond is denoted as 

We are considering the MG model as a representative example for the case that there is a gap in the spin excitation 
spectrum (and also in the charge channel). To give the result for the spectral density in the strong coupling limit 
J — > one has to find the modified quasiparticle residue Z(Q) in (J^). It can be simply derived from the correlation 
functions (see also Ref. |l2|) 

(si t ) = i[<#i|n,|#i> + <# 2 |n,|# 2 )] 



_> 



i\ n 1 / 1 x nJrl 



(fia n ) = (- 5 ) , <^Wi) = 2 ( "2 ) ' n ^ ' ( 43 ) 



in the following way 
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n=l 



3 1 + cos Q 
2 5 + 4cos2Q 



(44) 



The corrections for small J <C t may only be derived approximatively and we present two methods, projection method 
and variational procedure having different accuracy. 

A. Projection method 

First we calculate the spinon dispersion e s in the same way as it was done in the Heisenberg case in Sec. IV. 
But we should keep in mind that its applicability is less justified for the MG model than for the pure t-J model 
due to the much faster decay of spin correlation functions (compare with (|25|)). As before, we approximate 



(n 



(rii)(il,._ r /) which results in a constant energy shift from the term. Therefore, the first contribution to 



e s coming from J is merely determined by (fl' rr ,) (see Appendix B) resulting in 

£ s j(Q - 7r) = 2 J cos Q . 

We have a second contribution to e s from J' 

T 5 

es,/' {Q - ?r) = J' -4 cos Q + - + cos 2Q 

We see that for J' = J/2 the terms proportional cos Q cancel and we find 

"5 1 



e s (Q - 7r) = J 



cos 2Q 



(45) 



(46) 



(47) 



which is symmetric around tt/2. 

The spectral density is presented in Fig. 3. We see that in contrast to the t-J model the structures coming from 
Z(Q) (the holon branches) are much less pronounced, whereas square root singularities exist at the lower and upper 
edges of the spectrum. Their intensities are proportional to Z(k) or Z{k—Tr) at the lower and upper edges, respectively. 
Therefore, the square root singularity vanishes for k = ir at the lower edge. Furthermore, one can see that the low 
energy region for k between n/2 and n being empty in Fig. 2 is now filled with states. The spectrum becomes more 
symmetric around 7r/2 and the low-energy edge is given by the spinon dispersion e s (Q). The strong damping of the 
holon branch is due to the suppression of the singularity at the spinon Fermi edge (at Q = n/2 in Z(Q)). It is a 
universal feature for any ID magnetic state having a gap in the spin excitation spectrum. The suppression of holon 
weight was also found by Voitc3 for the Luther-Emery phase in the Luttinger liquid. The form of the spectral density 
in Fig. 3 resembles also roughly the exact diagonalization study in Ref. [12|. But a single bound state with a finite 
spectral weight that was obtained there, is missing in Fig. 3. That deficiency is due to the special projection procedure 
(|33|) which can only result in a continuous spectral density. Therefore, one has to go beyond the projection method. 



B. Variational ansatz 



Here we will use the set of string operators ( J18| ) as a set defining a variational wave function for the whole Hamil- 
tonian. Due to the knowledge of the exact ground-state ( (f2| ) all necessary matrix elements can be calculated without 
any further approximation. More precisely, we will diagonalize the Hamiltonian H = t + J + J' in the space spanned 
by the set of basis operators 



Vk,r 



1 

7l 



— oo 



(48) 



where L is the number of lattice sites and v m ,r was defined in (18). For that purpose one has to calculate the overlap 
matrix resulting in 



(49) 
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The kinetic energy part of the Hamilton matrix is given by: 

t-E^, = ({[v k ,r,t\,4,r>}) = \ (e ik (n r -r>-i) + e- ik (Q r _ r , . (50) 

The calculation of the exchange part of the Hamilton matrix is quite lengthy but straightforward. It shall not be 
given here in detail. To present the results we define a matrix Ef : 

^E* r ,=({[v k , r J+J>],vl r ,]) , (51) 
whose matrix elements are listed in Appendix C. The Hamilton matrix is then given by 

i=|i fe +pr> (52) 

and the matrix GF 

G r ,r> = {{Vk,r\ V k,r')) ( 53 ) 

can be found by solving the equation 

(u + iT + -Ejr 1 )£[ = g, T>0. (54) 

Finally the GF @ can be obtained by G(k,u>) = 2G , . 

The numerical results for J = and J = 0.4 at three different momenta are presented in Fig. 4. The curves for 
J = coincide with the analytic expression ( |28|]4~4| ). A number of 400 basis functions and a broadening of T = 0.05 are 
sufficient to reach the thermodynamic limit in contrast to the exact diagonalization method yielding only a sequence 
of 5-peaks. For J = 0.4 we can confirm the features found by the projection method, i.e. the low energy intensity 
between tt/2 and it, the symmetric spinon dispersion and the overdamped holon branch. In addition, the exchange 
terms produce two new features not present in Sec. VA: a resonance peak near zero energy and a bound state below 
the continuum. The resonance peak is visible near k = tt and becomes an antiresonance near k = 0. Careful inspection 
of the exact diagonalization dataE3 indicates also a very high peak at the resonance position for k — it and a small 
gap at k = 0, but a better understanding of the resonance/antiresonance feature is still required. 

The bound state is not visible in Fig. 4 due to the broadening T which is too large. Instead, we present in Fig. 5 the 
spectral weight of the lowest eigenstate w\ for k = it/ 2 and J = 0.4 in dependence on the number of basis functions. 
It is clearl*-|Seen that the weight tends to a constant value (wi w 0.1) in difference to the weight W3 of the third 
eigenstate.Ea At the same time, the separation e\ = £3 — E\ between the first and the third eigenvalues £1/3 stays 
finite for N — > 00 but the separation is very small (ei ~ 0.02 in units oli). For J = 0, both w\ and e\ tend to zero for 
N — > 00. That means that the bound state is connected with the presence of a gap in the spin excitation spectrum. 



VI. IDEAL PARAMAGNETIC STATE 



Such a state is realized for very high temperatures T, much larger than the exchange energy ksT 3> J. In that case 
spins at neighboring sites are completely uncorrelated. But the temperature is assumed to be lower than the Hubbard 
U such that the constraint of no double occupancy is preserved. Then the correlation functions become simply 



which results in 



Z(Q) = 



1 



~ I + cos Q 

The calculation of the spinon part (without frustration) gives 



e s (Q-7r) = 



J 



2 cos Q 



(55) 



(56) 







To calculate it one has to note that ( ]3q) i s no approximation in the present case. The effect of the (f2"r') terms can 
only be treated approximatively (see (p7|)) but it was checked by the variational method that its influence on the 
spectral function can be neglected. 

The information on Z(Q) and e s is sufficient to calculate the spectral function (Fig. 6). It is surprising that the 
strong singularities at the band edges survive despite the large temperature. The lower edge disperses according to 
the dispersion of the spinon (^6j) with a width proportional to J and has its minimum at k = n (in contrast to the 
frustrated case Fig. 3 with a minimum at k = 7r/2). But a peak connected with the holon dispersion proportional to 
t is not seen in Fig. 6. Such a peak appears in the finite temperature spectral function of the 2D t-J modeled and 
it can be expected since the first moment of the spectral function disperses according to tcosk. Its absence in ID is 
a nontrivial and unexpected result. It can be understood in the present context since the holon branch is strongly 
damped due to the suppression of the singularity in Z(Q) at Q — n/2. Apparently, that suppression is more strong 
in ( [Bq ) than in the frustrated case (Q) such that the holon branch is still visible in Fig. 3 but it disappears nearly 
in Fig. 6. One should note that the above result holds only in the region J <C k^T <C U. One may speculate that 
a further increase of the temperature such that the constraint of no double occupancy is lifted should lead to drastic 
changes in the spectral function. The strong singularities at the lower or upper band edges should disappear and a 
free dispersion should become visible. 

VII. CONCLUSION 

In conclusion we could derive analytic expressions for the spectral function of one hole in several magnetic states. 
The expressions are rigorous in the limit J — » 0, but our approach allows also to calculate the small J corrections. We 
analyzed the frustration and temperature effects. Results were given for the special frustration J' = J/2 with a gap 
in the spin excitation spectrum and for the ideal paramagnetic case. Both effects, frustration and temperature, lead 
to low-energy excitations between 7r/2 and 7r, and to a strong damping of the holon branches in the spectral function 
caused by the suppression of the singularity at the Fermi edge of spinons. The exchange terms in the MG model 
were found to be responsible for the finite weight of the lowest eigenstate and its finite, but small, energy separation 
from the rest of the spectrum, i.e. the bound state. The proposed scenario of holon branch damping seems to be a 
universal feature of frustration and temperature. Therefore, our results are of direct importance for photoemission 
experiments on strongly frustrated ID compounds like CuGe03, for instance. However, edge-shared-cuprate chains 
have a smaller energy scale and less ideal ID behavior in comparison with corner-shared compounds ,t3 which hinders 
direct comparison with experiment. But it cannot be excluded that a small frustration is also present in SrCu02 such 
that our study gives one possible reason, whs: no real, separate holon branch could be observed in the experimental 
spectra of SrCu02 between k = and 7t/2.lI In the spin gap case we found a very small energy separation of the 
bound state from the continuum such that it is nearly impossible to detect it in a photoemission experiment. 
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Appendix A: Spinon dispersion of the Heisenberg case 



In this Appendix we outline the main steps to derive the spinon dispersion of the pure t-J model using the projection 
method. For long chains of X-operators it is convenient to introduce the notationscJ 



\ yaai v a l a 2 

2^ A »l A n 2 A n 3 



(ni\n 2 \n 3 \ ■ ■ ■ |n r _i|n r 



and 



X™ 1X % a2X % a3 ■ ■ ■ KZl^Kr = (ni\n 2 \n 3 \ . . . \n r ^\n r ) 



which means especially that ( p"8| ) may be rewritten as w m , r = (m\m + g\m + 2g\ . . . \m + r]. In such a notation we 
obtain for the commutation with the Heisenberg Hamiltonian 



= ~i E X ™+9 X ™ = -i E(™ + 9\m] , 



(A.l) 



x af3 f 



E (x%X%_ g - Xl\ g X^) = J - £{(to|to + g\ - (m + g\m\} , 



(A.2) 



and then 



— {(m\m — g\m + g\m + 2g| . . . \m + r] — (m — (?|to|to + g|m + 2<?| . . . \m + r] 

— (to + g\m + 2g\ . . .\m + r] — (m| . . . |m + r — g\m + r + g\m + r]} , 



(A.3) 



where g = sign(r). In deriving (A.3) it is important that the commutation of the "inner" operators X^f +l with I < r 
do not give rise to additional terms since the corresponding sums cancel each other. That is a direct consequence of 
one-dimensionality. 

Now, we consider the holon contribution to e s coming from v'{, 

({< 8 .4^ 9 '}}= 2 ^(fc-fcOE ei9 ' r '^ r+ifc(r " r ' ) ^) ( A - 4 ) 

r,r' 

with 

£l r r i — (rn\ . . . \m + r — g\m + r + g\m + r\m + r — g'\ . . . \m + r — r ) , (A. 5) 



and g — sign(r), g' — sign(r'). We see that in general (£l" r >) depends both on r — r' and on r. Eqn. (A^j) can also 
be written as 



{vl q ,vi, q ,}) = 2n5(k-k>)J2e lik - q)l Si 



(A.6) 



with 



+00 



a= E 



-i(q-q')r /q/i 



r.r — u 



r— — 00 



Due to the slow decay of spin correlation functions in the ID Heisenberg state, one can expect that the main 
contr ibution to Si comes from regions where |r| ^> \l\. There holds g = g' and we may rewrite and approximate 
©by 

(n^, r _ I ) = ((o|...|0(r+ 5 |r))«<n / )(n 1 ). (aj) 

Then, the explicit dependence on r drops out and we obtain 
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i.e. a simple constant shift of the energy e s . 

The contribution of the v' k term to the spinon dispersion is determined by the sequence of spin operators 

fl' r r i = (m|m — g\m + g\ . . . \m + r|m + r — g \ . . . \m + r — r ) 



(A. 



(A.9) 



instead of (A. 5). The expectation value of that term has to be calculated for the magnetic system without holes. In 
difference to (f2" r ,), it depends only on r — r' without further approximation 



(n^. Ijr ) = <n{ i0 ) = <(o|-i|i|...|0) ) (i>o), 

and o) = o)- For I = 0, 1 it can be expressed through pair correlation functions 



(A.10) 



For large I > we may expect 



Ko> = J + 2(S S 1 ) + (S S 2 > 



<n! ) « (n I+1 ) 



(All) 



Using this approximation we obtain the following contribution to the spinon dispersion e s which stems from the v' k 
term 



\ ({< q > v l, q >}) = 8 ^ l') 5 ( k fc 'X(fc - Q)Z(k -q + n) 



with 



J 



+ 00 



e' s (k)Z(k + ir) = - £ e ~^,o 



(A.12) 



(A.13) 



After some algebra we find 

Ae' s (Q-7r)Z(Q)/J= {n' 0>0 ) -2<fii) -2cosQ 

where 



sin QY(Q) , 



Y{Q) = £ 2(-l)' Bin(Q0<nj> = ^ / <fcZ(«) 



z=i 



Q — k Q + k 

COt h COt 



The complete expression for the spinon dispersion follows from 

e s (Q — 7r) = J cos Q + e^(Q — n) + const 
and is given in Eqn. (B9|) neglecting the constant energy shift. 



(A.14) 



(A.15) 
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Appendix B: Majumdar-Ghosh model with projection method 



The commutation with the frustration Hamiltonian (|14|) is very similar to ( |A.3[) . It gives 

J' 



Urn. . r i J 



1 m,r ' rn,r J 



(B.l) 



where 



'm r = {(w|m — 2g\m + g\ . . . \m + r] — (m — 2g\m\m + g\ . . . \m + r] 

— (m\m + g)(m + 2g\ . . . |m + r] + (m\m + g\m — g\m + 2g\ . . . \m + r] 

— (m|m — g\m + g\ . . . \m + r]} , 



,(4) =/_ 



= { — (m\ . . . \m + r — g\m + r + 2g\m + r] 

— (m| . . . \m + r — 2g|rn + r + <?|m + r — <?|m + r] 
+ (m\ . . . \m + r — g\m + r + g|m + r]} . 



(B.2) 



Again we see that all terms coming from commutations at "inner" operators cancel. But now the motion of spinons 
becomes more complicated. The same considerations that show the absence of dispersion from the v'^ r term are 

applicable to Vm]r- 



Taking into account that 



n > , 



we obtain the following contributions to e s = e s j + e s j/ (we drop dispersionless terms). From J comes 

J 



e sJ (Q-7T)Z(Q) 
Z'(Q) 



Thus 



[2cosQZ(Q) + Z'(Q)] , 

z, 

3(-±+cosQ) 11 

5 ^ + - + - = Z(Q 2cosQ 

5 + 4cos2Q 4 8 W) * 



£sj(Q — 7r) = 2 J cos ( 



has the same form (JO]) that we have assumed for the t-J model. 

For the term that comes from the left distorted end of v m r due to J' we have 



\\ = I /o (3) \ 

^"m,ri' KJ m+r 1 —r2,r 2 J / — 2 \ r i~ r 2 / ' 



(B.3) 



(B.4) 



n 



(3)\ _ 
2n / 



(3) 



-4( -- ) , n>l, 



(3) 



The contribution from J' is 



e sJ -(0-7r)Z(g) 



n 



(3) 



-4Z'(Q) + ^ + ^cosQ 
4 



z(Q) 



j' 



icosQ + 2(- + cos2Q) 



(B.5) 



and we obtain 
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Appendix C: Matrix elements of the variational basis set 

The matrixelements of E* , in the neighborhood of r, r' = are given by: 







-3 


-2 


-1 





1 


2 


3 


r 


















-3 




3/8 





3/16 


-3/32 


-9/32 


3/32 


3/16 


-2 







3/8 


-9/16 


3/8 





-3/8 


3/32 


-1 




3/16 


-9/16 


3/8 


-3/16 


3/8 





-9/32 







-3/32 


3/8 


-3/16 





-3/16 


3/8 


-3/32 


1 




-9/32 





3/8 


-3/16 


3/8 


-9/16 


3/16 


2 




3/32 


-3/8 





3/8 


-9/16 


3/8 





3 




3/16 


3/32 


-9/32 


-3/32 


3/16 





3/8 



One has two different regions in the matrix. The first one is defined for r > 0, r' > and r > r' + 2 where we have 
the matrix elements: 



r = 2n, r' = 2m or r = 2n + 1, r' = 2m + 1 : E'f. r , = 



r = 2n + 1, r' = 2m : 
r = 2n, r' = 2m — 1 : 

and the second one for r' < — 2, r > 2 with 



E x 



3 ( l\n-m 

2) ' 

I \ n— m+1 



8 ( 2) 



(C.l) 



E* r , - --(Cl r - r ') . 



(C.2) 



There are special matrix elements along the diagonal (3/8) and along the side diagonal (alternatively -9/16 or 0) and 
also for the two lines (n > 1): 



The matrix is filled by 



E. 



2n,0 



3 

16 



E2n,-1 — ^2n+l,-l — Yg 



J^X J^X J^X 



(C.3) 
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Figures 

Fig. 1: Comparison of the spinon dispersion e s (Q — n) as calculated from the projection method (broken line, 
J = 0.4) with 2 J cos Q (full line). 

Fig. 2: Spectral density of the t-J model for J = 0.4 and t = 1. 

Fig. 3: Spectral density of the frustrated t-J model (J = 0.4 and t = 1) at the special frustration J' = 0.5 J (using 
the Majumdar- Ghosh wave function) within the projection method. 

Fig. 4: Spectral density of the Majumdar-Ghosh model A(k,uj) for three different momenta k/ir and t = 1, J = 0.4 
(full lines) or J = (dashed lines) with a variational set of 400 basis functions and a broadening of T = 0.05. 

Fig. 5: Weight w\ and energy separation ei = E3 — E\ of the lowest eigenvalue E\ at k/ir = 0.5, t = 1, J = 0.4 
(full lines) as a function of the inverse number of basis functions l/N. The dashed lines are the weights W3 and the 
energy separation = E§ — E3 of the third eigenvalue £3 . 

Fig. 6: Spectral density of the t-J model (J = 0.4 and t — 1) in the ideal paramagnetic state. 
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